SALAMANDRA 50(3) 


ISSN 0036-3375 


147-154 30 October 2014 


Assessing factors involved in determining fluctuating asymmetry 
in four insular populations of the Balearic lizard Podarcis lilfordi 

Mario Garrido & Valentin Perez-Mellado 

Department of Animal Biology, University of Salamanca, Campus Unamuno, 37007, Salamanca, Spain 

Corresponding author: Mario Garrido, e-mail: gaiarrido@usal.es 

Manuscript received: 10 December 2013 
Accepted: 23 April 2014 by Philipp Wagner 


Abstract. In animals, developmental stability is frequently assessed by the level of fluctuating asymmetry. Several environ¬ 
mental and genetic stress factors can increase the developmental instability in a population. Anyhow, the use of fluctuat¬ 
ing asymmetry as a measurement of developmental instability and its relationship to other measurements of genetic and 
environmental stress, remain controversial. We studied this subject in the Balearic lizard, Podarcis lilfordi, and examined 
four populations inhabiting different coastal islets of Minorca (Balearic Islands, Spain). These populations show clear eco¬ 
logical and genetic differences that allowed us to study factors potentially determining the level of fluctuating asymmetry. 
We examined the following predictions: (i) fluctuating asymmetry will increase on smaller islands; ( 2 ) fluctuating asym¬ 
metry will increase on islands with a higher lizard density; ( 3 ) fluctuating asymmetry will increase on islands with more 
heavily parasitized lizards; and ( 4 ) populations with higher inbreeding levels would be more susceptible to environmental 
stress. Our results would partially support the fourth prediction of a higher fluctuating asymmetry in populations with 
lower genetic variability and, consequently, a higher inbreeding level. Probably, other environmental factors modulate this 
relationship. 
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Introduction 

Developmental stability is viewed as an individual’s abil¬ 
ity to buffer its development against random perturbations 
(van Dongen & Lens 2000a). It is frequently assessed 
by the level of fluctuating asymmetry (FA), which corre¬ 
sponds to the small and random differences between the 
right and left sides of an otherwise bilaterally symmetri¬ 
cal character (van Valen 1962). Several environmental and 
genetic types of stress can increase the degree of develop¬ 
mental instability. Fluctuating asymmetry today tends to 
be considered a population parameter, playing no role at 
individual level (Palmer 1999). However, it is also used as 
an indicator of an individual’s quality (M0ller & Swad¬ 
dle 1997, VAN Dongen & Lens 2000b). 

This study tries to approach the subject of FA using the 
Balearic lizard, Podarcis lilfordi, as a model species. We ex¬ 
amined four populations of P. lilfordi inhabiting different 
coastal islets around the island of Minorca (Balearic Is¬ 
lands, Spain). Although ectotherms are generally consid¬ 
ered to be good indicators of environmental stress (Leary 
& Allendorf 1989), there are few studies in the develop¬ 
mental stability of reptiles. In addition, lizards have many 
metric (such as limb or head dimensions) and meristic 
(such as scales, femoral pores or subdigital lamellae) paired 
characteristics that are relatively easy to measure (Soule 
1967), and they are excellent models for ecological or eco- 


morphological studies (Pianka 1986, Losos 2009). At 
population level, the small islets around Minorca present 
ideal scenarios to study the effects of genetic and ecological 
factors that may influence development stability, because 
they host isolated lizard populations in relatively simple 
ecosystems. Also, developmental instability measured by 
FA is, presumably, more evident in small island popula¬ 
tions because they may be exposed to more environmental 
(Mac Arthur & Wilson 1967) and genetic stress (Young 
et al. 1996, Frankham 1998, Crawford et al. 2001). 

In the absence of human intervention, gene flow be¬ 
tween islands is limited or absent. Crnobrnja-Isailovic 
et al. (2005) argued that isolated populations were more 
prone to be impacted by conditions such as bottleneck ef¬ 
fect, long-term strong selection, low effective population 
size, or crowding effects, which may lead to higher in- 
breeding levels, periods of intense directional selection, 
and extreme homozygosity. In addition, Moller & Swad¬ 
dle (1997) argued that all these factors “have been shown 
to break down the balance of a genome and reduce the 
ability of the organism to buffer its developmental path¬ 
ways against the production of random errors”. Hence, the¬ 
ory suggests a positive correlation between FA and levels 
of inbreeding as it would increase homozygosis (Moller 
& Swaddle 1997). 

In addition, small populations may also withstand more 
environmental stress through, for example, competition for 
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limited resources and/or co-existence with predator spe¬ 
cies (Mac Arthur & Wilson 1967, Crnobrnja-Isailovic 
et al. 2005). Due to increased genetic stress, small popu¬ 
lations may be more susceptible to environmental stress. 
While not-inbred populations do not suffer from certain 
changes in the environment, inbred populations perceive 
them as stressful, which indicates a causal link between ge¬ 
netic and environmental stress (Kristensen et al. 2006). 
Such connection can lead to synergetic effects between 
both types of stress, resulting in environmental conditions 
that would normally be benign becoming detrimental 
(Kristensen et al. 2006). To assess the impact of compet¬ 
itor or predator species, we deliberately selected popula¬ 
tions where we have four combinations of presence or ab¬ 
sence of rats and gulls available for study. These two species 
are often perceived as significant agents in the alteration of 
many ecological communities, particularly on islands (Pe¬ 
rez-Mellado 2002, OrO & MARTfNEZ-AsRAfN 2OO7). 
However, Mus spretus Lataste, 1883 is present on Colom 
Island, and even though mice are not normally recognized 
as predators of Mediterranean lizards, they are known as 
predators of Teira dugesii (Milne-Edwards, 1829) on At¬ 
lantic islands (P. Wagner, pers. comm., 19 April 2014). 

Among non-genetic factors, parasite load and popula¬ 
tion density have been the most widely studied factors in 
relation to FA. Moller & Swaddle (1997) and Moller 
& Manning (2003) presented large sets of evidence from 
a range of vertebrate taxa and concluded that population 
density could result in higher levels of FA, presumably 
generated by the reduced availability of resources at high¬ 
er densities and increased competition. However, other 
authors did not find such a relationship between popula¬ 
tion density and FA (Wiener & Rago 1987) or in general 
(Kruuk et al. 2003). 

Another potential source of environmental stress is 
parasitism, which is considered to be detrimental to host 
individuals. It also has been linked with the occurrence of 
asymmetry in various traits (Mgller 1996 and referenc¬ 
es therein). Nevertheless, little is as yet known about the 
possible causative nature of the relationship between par¬ 
asite load and developmental instability (Polak 1997). In 
two different reviews, Mgller (1996, 2006) supported the 
hypothesis that there was a positive relationship between 
parasite load and FA. He also warned that the link between 
both variables were expressed in various ways, did not im¬ 
ply causation, and that the damage caused to the host was 
related to the parasites virulence. As Brown & Brown 
(2002) proposed, developmental instability itself would 
affect the immune system and cause greater susceptibility 
to parasites. They also proposed that poor environmental 
conditions (e.g., poor nutrition, pollutants) would promote 
developmental instability, and it follows that these might 
also promote the parasite load of hosts by increasing their 
exposure to parasites, perhaps as a result of positive effects 
on local parasite abundance. Brown & Brown (2002) con¬ 
cluded that these possibilities have not been ruled out in 
most of the studies that have shown a relationship to ex¬ 
ist between parasite load and asymmetry (Moller 1996, 


Moller & Swaddle 1997). In theory, populations that are 
exposed to more stressful conditions should exhibit high¬ 
er levels of parasite load and FA. Notwithstanding, sever¬ 
al studies did not find such a correlation (Bergstrom & 
Reimchen 2005 and references therein). 

In lizards, the isolation of insular populations has been 
shown to lead to increased inbreeding and genetic drift 
(Madsen et al. 1996, Madsen et al. 2000), and the size of 
an island has been associated with the degree of FA occur¬ 
ring there (Soule 1967, Sarre & Dearn 1991, Sarre 1996, 
Vervust et al. 2008, Bancila et al. 2010). However, no 
such increase in FA was found in small island populations 
of Podarcis muralis (Crnobrnja-Isailovic et al. 2005). 
At individual level, some studies related the asymmetry of 
hatchlings to temperatures during egg incubation (Qualls 
& Andrews 1999, Ji et al. 2002, Zhdanova & Zakha¬ 
rov 2006), physiological performance (MartIn & Lopez 
2001, Lopez & MartIn 2002, Vervust 2008, but see also 
Warner & Shine 2006), or mating and reproductive suc¬ 
cess (MartIn & Lopez 2000, 2006, but see Lailvaux & 
Irschick 2006). 

The objective of our work is to study the influence of 
genetic and non-genetic (density, competitor species and 
parasitism) factors on the degree of FA in both individuals 
and populations. We predict that FA will increase: (1) on 
smaller islands; (2) on islands with a higher lizard density; 
and (3) on islands with more heavily parasitized lizards. We 
also expect that populations with higher inbreeding levels 
would be more susceptible to environmental stress, that is, 
the impact of competitor species should increase with in- 
breeding levels. 

Material and methods 

Species under study 

The Balearic lizard, Podarcis lilfordi (Gunther, 1874) 
(Squamata, Lacertidae), is a medium-sized lizard endemic 
to the coastal islands and islets of Mallorca, Minorca and 
the Cabrera archipelago (Balearic Islands). Lizards were 
collected by means of noosing. For each individual, we re¬ 
corded its age as adult or juvenile and sex, SVL, tail length 
(if original), and body weight. Because of the low number 
of juveniles captured, only data from adult individuals 
were analysed. Body condition was evaluated from the re¬ 
sidual of the regression of log-transformed body weight 
versus log-transformed SVL. 

We quantified the prevalence and intensity of infection 
by two different parasites, mites and haemogregarines. 
Both kinds of parasites were previously studied (Garrido 
& Perez-Mellado 2013a, b). Mites of P. lilfordi were pro¬ 
visionally identified as belonging to the genus Ophionys- 
sus, while blood parasites were regarded as representing 
the haemogregarines group (Garrido & Perez-Mellado 
2013a, b). Fctoparasites were counted at the site of capture 
with a 5X eyeglass, inspecting the whole body surface of 
the lizards. To assess the extent of blood parasitism we ob¬ 
tained blood samples by making a slight longitudinal cut 
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Table 1. Genetic data and environmental factors of each population under study: NFS: total number of polymorphic sites; Pi: nucle¬ 
otide diversity; k: average number of nucleotide differences; S: island area (in hectares); Rats: presence (-I-) or absence (-) of the ship 
rat, Rattus rattus-. Seagulls: presence (-I-) or absence (-) of seagull breeding colonies; D: Biotic capacity index; B: niche breadth index 
of vegetation and substrate cover; lizard density (mean ± SE lizards ha“'; Biome: lizard biomass (kg ha'). Source: Perez-Mellado 
et al. (2008). 


Population 

NFS 

Pi 

k 

s 

Rats 

Seagulls 

D 

B 

density±SE 

Biome 

Aire 

1 

0.000200 

0.50 

29.80 

- 

+ 

6.10 

5.58 

4098.60±586.60 

32.54 

Colom 

4 

0.00095 

2.07 

51.14 

+ 

+ 

7.62 

9.38 

1615.50±294.93 

11.02 

Binicodrell 

0 

0.00 

0.0 

0.29 

- 

- 

0.57 

3.84 

1180.50±264.45 

7.59 

Tosqueta 

0 

0.00 

0 

0.29 

+ 

- 

0.57 

4.91 

1703.00±501.41 

11.29 


in the dorsal proximal side of the tail with a sterile scalpel. 
The blood of drop emerging allowed us to obtain a blood 
smear in situ. Specimens were always released at the site of 
capture. Blood smears on microscopic slides were air-dried 
in the field. In the lab, these slides were fixed with abso¬ 
lute methanol for lo minutes and then stained with modi¬ 
fied Giemsa for 20 minutes. Samples were analysed using 
an optical microscope at 400X. The only blood parasites 
identified were haemogregarines. The density of haemo- 
gregarines was estimated on the basis of a total of 2,000 
counted cells per sample. Their prevalence was estimated 
as the percentage of infected individuals on each island. 
Within each year, individuals were made recognisable 
marking them with colour pens before their release. Liz¬ 
ards already examined in previous years could be recog¬ 
nized by the scar of the cut in the tad. Balearic lizards were 
studied on four different islands: Aire (n = 133 females, 
183 males), Colom (n = 89 females, 121 males), Binicodrell 
(n = 19 females, 31 males), and Tosqueta (n = 8 females, 
33 males). In all cases, individuals were sampled during the 
springs and summers of 2007, 2008, 2009 and 2010. The 
relatively simple ecosystems of these islands allowed us to 
identify the ecological variables that differed between pop¬ 
ulations. We deliberately selected these four populations, 
because the ecological and genetic differences between 
them were already known (Perez-Mellado et al. 2008, 
see Tab. 1), allowing us to study the factors potentially de¬ 
termining their levels of FA. 

FA assessment 

We measured one metric (hind limb length, HLL) and six 
meristic traits from the left and right sides: the number of 
femoral pores on the thigh (FEM), the number of subdig¬ 
ital lamellae under the fourth toe (LAM), the number of 
supraocular scales (SPO), the number of supralabial scales 
(SPL), the number of infralabial scales (INL), and the 
number of submaxillar scales (SMX). To verify whether the 
traits met the requirements for assessing FA (see below), 
we took three repeat measurements of them in a sample of 
30 lizards. Hind limb length was measured from the inser¬ 
tion point of each leg into the trunk to the base of the claw 
of the fourth toe while keeping the leg fully extended. We 


took the measurements with digital callipers to the nearest 
0.01 mm. Meristic traits were counted with a 5X eyeglass at 
the site of capture. 

For metric and meristic traits in which the absolute dif¬ 
ferences of right minus left (|R - L|) sides are large (4-5 
units). Palmer (1994) and Palmer & Strobeck (1986) de¬ 
scribed a test to establish whether the between-sides vari¬ 
ation was significantly larger than the measurement error 
(ME). The test is a two-way ANOVA (sides x individuals) 
and should be conducted routinely as part of any study of 
FA (Palmer & Strobeck 1986). This procedure tests for 
the significance of all between-sides variation relative to 
measurement error, including antisymmetry (AS) and FA. 
One advantage of factorial ANOVA procedures is that the 
significance of directional asymmetry (DA) can be tested 
at the same time as that of FA relative to measurement er¬ 
ror. For meristic traits, where the difference between sides 
is small (1-2 units) and where a low level of scoring er¬ 
ror may exist, ANOVA techniques do not work properly, 
though. Palmer (1994 and pers. comm.) therefore recom¬ 
mended applying a likelihood ratio or G-test (Sokal & 
Rohlf 1995). If the test produces a statistically significant 
result, it still needs to be confirmed that the mean differ¬ 
ences in the character state between both sides (|(R - L)/n|) 
is greater than the mean difference between all measure¬ 
ments (|Mx - My|/n), where ‘n’ is the number of replicate 
counts. In these cases, DA must be estimated with a one- 
sample t-test of the mean (R - L), because departures from 
normality can mask DA (Palmer 1994). 

Once tested for DA and ME, tests for the presence of AS 
should also be conducted in studies of EA. Typically, AS 
shows a platykurtic or bimodal frequency distribution in 
a given population, while DA presents skewness (M0ller 
& Swaddle 1997). To check for kurtosis and skewness, we 
applied the Anscombe-Glynn and DAgostino tests, re¬ 
spectively (Zar 2010). Finally, Palmer & Strobeck (1986) 
suggested looking for a possible dependence of FA to body 
size, correcting for this dependence if necessary. Palmer 
(1994) pointed out that body size might reflect some as¬ 
pect of body condition. Thus, checking for trait size may 
partially control the extent of influence by body condition 
and yield asymmetry estimates that are condition-inde¬ 
pendent. Accordingly, we chose to do without a correction 
for size dependence, because body size could be correlat- 
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Table 2. Results of G-tests and the protocol (see Palmer 1994) to assess FA relative to ME from meristic traits whose average differ¬ 
ence between body sides is low. Last row: results for one-sample t-test, testing for departure from the expected mean of (R - L) = zero. 
This method allowed discarding the presence of DA in the different characters. 



FEM 

LAM 

SPO 

SPL 

INL 

SMX 

x±SE of |R- L| 

0.75±0.04 

0.83±0.05 

0.07±0.01 

0.12±0.02 

0.19±0.02 

0.03±0.01 


G, = 34.89 

G, = 6.54 


G , = 7.72 

G, = 4.95 


Palmer’s G-test 








P < 0.001 

P < 0.001 


P < 0.05 

P < 0.05 


x|Ml - M2| 

0.03 

0.49 


0.01 

0 


x|(R-L)/2| 

0.074 

1 


0.14 

0.07 


FA > ME 

yes 

yes 

yes 

yes 

yes 

yes 


L„ = -0.35 

L,. = -1-01 

t = - 0.93 

L,. = -0.12 

L„ = 0.19 

t„, = 0.89 

Fisher’s t-test 








P = 0.73 

P = 0.31 

p = 0.35 

P = 0.91 

P = 0.83 

P = 0.37 


ed with factors under study and such a correction might 
mask their association with FA (Palmer & Strobeck 1986, 
Bancila et al. 2010). 

For our purpose, we used a composite FA index (CFA), 
as the analysis of composite indices is more powerful than 
analyses based on single traits (M0ller 2006, see discus¬ 
sion in Leung et al. 2000). In addition, analyses that com¬ 
bine information across traits should prove to be more 
reliable detectors of stress (Leary & Allendorf 1989, 
Watson & Thornhill 1994). Therefore, after all traits had 
been analysed for FA, they were combined per individu¬ 
al in the composite CFA2 following the recommendations 
BY Leung et al. (2000). Briefly, the FA value of each trait 
was divided by its average |FA|, so that all traits contribut¬ 
ed equally to the CFA, and then |FA| values were summed 
across traits for each individual so that each individual 
achieved a composite FA score. 

Statistical analysis 

First, we ran a GLM model with CFA2 as the dependent 
variable, and island and mite and haemogregarine loads as 
independent variables. To avoid spurious differences that 
are actually caused by sex-dependent differences, we in¬ 
troduced this factor to the model. CFA2 data were Box- 
Cox transformed to obtain homoscedasticity. Then, as we 
obtained differences between insular populations, we ran 
another GLM model, introducing environmental (rats, 
seagulls and density) as well as genetic factors of popula¬ 
tions (NPS, Pi and k; Tab. 1). Thus, we looked for the fac¬ 
tors leading to differences in FA between islands. We tested 
for co-linearity of the explanatory variables and discarded 
any variable with a variance inflation factor (VIF) of > 3 
(ZuuR et al. 2010). Genetic variables showed VIF values 
higher than this threshold, for which reason Pi and k were 
deleted and NPS was introduced to the model, producing a 
final set of four variables. 

The statistical analyses were carried out using the soft¬ 
ware R environment (ver. 2.12.1, R Development Gore 
Team, 2010). We started out with saturated models, con¬ 


sidering effects of all independent variables as well as all 
interactions between them. To determine the minimum 
adequate model (Grawley 2007), we used a backward 
stepwise regression on the basis of deletion tests (Craw¬ 
ley 2007). Non-significant factors and interactions were 
removed from the models until we obtained a minimal ad¬ 
equate model, taking AIC (Akaike’s Information Criterion) 
into account at each step. Post-hoc comparisons were made 
using the multcomp package (Hothorn et al. 2009). Effect 
sizes reported here are partially eta-squared (qp that can 
be interpreted similarly as R^ (Cohen 1988). Like Rh qj is a 
dimensionless measure of effect size that is independent of 
the degrees of freedom used in the analyses (Tabachnick 
& Fidell 2000). 

Results 

Assessing FA 

All meristic traits show an average difference between sides 
of < 1 (Tab. 2). We therefore applied Palmer’s recommended 
G-test to see if ME was larger than FA (Tab. 2). The result 
of the G-test was significant in that the mean of |(R - L )/2| 
was greater than the mean of |Mi - Maj for all traits except 
SPO and SMX where no G-test could be applied as there 
was no variation in the repeated measurements (100% re¬ 
peatability of measurements; Tab. 2). For these traits, we 
obviously concluded that variation between sides could 
not be due to ME. Moreover, the normality of all meristic 
traits was tested to be exclusive of DA (see above), and we 
found that distributions were normal for all meristic traits 
(Tab. 2). For HLL (metric trait), the two-way ANOVA rec¬ 
ommended by Palmer showed that DA was not present 
and that the variation due to non-directional asymmetry 
was higher than the variation due to ME (Tab. 3). 

A visual inspection of data plots revealed no evidence 
for a clustering of vectors of shape asymmetry that could 
have suggested the presence of AS. Nevertheless, signifi¬ 
cant levels of skewness for the distribution of right minus 
left-differences were found for SMX, but not for the re¬ 
maining traits (Tab. 4). All traits showed a distribution that 
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Table 3. Two-way mixed model ANOVA: testing the significance of between-side individual interaction (non-DA), between-side 
variation (DA), and between-replicates variation (ME) in Podarcis lilfordi. The asterisk (*) indicates significant factors; NS indicates 
no significant factors. 




Df 

Sum Sq 

Mean Sq 

F 

P 


side 


1 

0.00 

0.001 

0.0085 

0.9265 


Individual 


25 

1090.96 

43.638 

274.8742 

< 0.001* 


SideTndividua] 


25 

28.37 

1.135 

7.1472 

< 0.001* 


Residuals 


104 

16.51 

0.159 




Non-DA (FA and AS) 

DA 

Between-replicate 

Trait 

MS 

df 

F 

MS 

df 

F 

variation (MS^) 

HLL 

1.135 

25,104 

7.14* 

0.001 

1,25 

O.OOOO*^® 

0.159 


was significantly leptokurtic (Tab. 4). Since there were no 
platykurtosis or bimodality, we eliminate the possibility of 
AS acting as a confounding factor in our tests for FA. 

To sum our results up, a significant degree of FA was 
detected in all characters examined in P. lilfordi, except for 
the skewness found in SMX, probably due to the test being 
limited to just two individuals that were asymmetrical for 
this trait and as a consequence of the small extent of vari¬ 
ation in this trait. We decided to discard SMX in further 
analyses. 

FA analyses 

In the GLM model with parasite loads, sex and island, all 
interaction terms were insignificant (all P > 0.05), so they 
were excluded from the model. Moreover, we did not de¬ 
tect differences in FA between males (mean ± SE = 6.01 
± 0.34; n = 249) and females (mean ± SE = 4.68 ± 0.28; 
n = 368; F^ = 1.90, P = 0.17), and the effect size was very 
small (qj = 0.004). Also, no correlation was found for the 
haemogregarine load (F^ - 0.03, P = 0.87). Both were ex¬ 

cluded from the final model. 

Mite loads showed no significant correlation with CFA2 
(F^ - 2.52, P = 0.11) either, but were retained in the mini¬ 

mum model. We found strong differences among islands 
(F3 ^3^ = 19.34, P < 0.001), but the effect size was moder¬ 
ate (q^ = 0.094). Post hoc tests revealed that Tosqueta 
(mean ± SE = 12.68 ± 1.07, range = 23.89) showed statisti¬ 
cally higher rates of asymmetry than Binicodrell (mean ± 
SE = 8.15 ± 1.47, range = 60.52; t = 4.05, P < 0.001), Colom 
(mean ± SE = 5.91 ± 0.38, range = 27.89; t = 6.63, P < 0.001), 
and Aire (mean ± SE = 5.14 ± 0.26, range = 25.10; t = 7.52, 
P < 0.001). Binicodrell individuals were found to be dif¬ 
ferent from those living on Aire (t = 2.68, P = 0.04). Bini¬ 
codrell did not show differences from Colom (t = -1.80, P = 
0.26). Lizards from Colom were not more asymmetric than 
those from Aire (t = 1.31, P = 0.54). 

The GLM model with environmental and genetic fac¬ 
tors exhibited no significance for either the interactions or 
the presence or density of seagulls (all P > 0.05). Only NPS 
(El = 46.27, P < 0.001) and the presence of rats (F^ = 

10.94, P - 0.001) were related to CFA2. Thus, populations 


Table 4. Skewness (DAgostino test) and kurtosis (Anscombe- 
Glynn test) estimates for each measured trait. The asterisk (*) 
indicates significant factors. 


Trait 

Skewness 

P 

Kurtosis 

P 

LMP 

0.69 

0.48 

3.81 

< 0.001* 

FEM 

0.1813 

0.26 

4.080 

< 0.001* 

LAM 

-0.0022 

0.99 

4.099 

< 0.001* 

SPO 

-2.31 

0.02 

11.54 

< 0.001* 

SPL 

0.19 

0.24 

9.97 

< 0.001* 

INL 

1.59 

0.11 

7.37 

< 0.001* 

SMX 

2.44 

< 0.001* 

14.69 

< 0.001* 


living without the company of rats and with more genetic 
variability were less asymmetric. 

Discussion 

We used CFA2 to test some hypotheses on FA. We did not 
find any differences between sexes in the levels of FA in 
contrast to some previous studies where males were more 
asymmetric than females (M0ller & Swaddle 1997). 
However, our result might be a consequence of the chosen 
sampling design. According to the theoretical assumptions 
of FA, we can conclude that both sexes are exposed to simi¬ 
lar levels of stress. 

FA levels were significantly different between all four is¬ 
land populations, with the exception of Aire versus Colom 
where similar degrees of asymmetry were found. Although 
we had expected to find higher FA levels on smaller islands 
due to increased levels of inbreeding, there was a lack of 
difference between Binicodrell and Colom, probably due 
to unequal sample sizes and a wide range of FA values on 
Binicodrell. Density has been pointed out as severely af¬ 
fecting developmental stability and as a consequence in¬ 
creasing FA (see, for example, M0ller & Swaddle 1997, 
M0LLER & Manning 2003), and very few studies have 
questioned this relationship (e.g., Kruuk et al. 2003). In 
our study, this relationship is not as obvious as the theory 
would suggest. GLM showed no influence of the density 
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on CFAi, but this could be due to our small sample size of 
just four island populations. Aire Island had a much high¬ 
er population density, yet produced the lowest FA values 
along with Colom. Among the other islands, with similar 
population densities, FA differed substantially (Tab. i). 

Immune system and parasitism also could be linked to 
FA in poor environmental conditions (M0ller & Swad¬ 
dle 1997 and see Introduction). But in our case, neither 
mites nor haemogregarines seem to be good predictors of 
the degree of asymmetry. Parasitism not only depends on 
immune system efficiency, but also on other ecological and 
environmental factors (Biaggini et al. 2009). It even de¬ 
pends on the behaviour of individuals within a population 
(Garrido & Perez-Mellado 2013a, 2013b, and references 
therein). Binicodrell individuals provide a good example of 
how, even in the absence of ectoparasites, we can find high 
FA levels. Also, negative results could be a consequence of 
our sample design. Here we studied two kinds of parasites 
with different degrees of virulence (see Material and meth¬ 
ods) and just four populations with very different intensi¬ 
ties of infection. 

According to our results obtained from four insular 
populations, genetic stress would be a more reliable in¬ 
dicator of the degree of asymmetry. Theory predicts that 
small populations may suffer more genetic stress than large 
populations due to the erosion of genetic variation (Young 
et al. 1996, Frankham 1998, Crawford et al. 2ooi).This 
loss of genetic variation and deeper inbreeding depres¬ 
sion can cause declining fitness (Crnokrak & Roff 1999, 
Frankham 2005, Allendorf & Luikart 2007). Genetic 
variability was significantly higher in larger populations 
(Tab. 1). Amongst the two less asymmetrical populations, 
Colom and Aire, no differences in CFA2 were found, but 
Colom showed a comparatively higher genetic variability 
(Perez-Mellado et al. 2008). Despite this exception, our 
results are consistent with other studies and prove the influ¬ 
ence of genetic variability on FA levels. We propose that ge¬ 
netic variability affects stability when it has declined below 
a certain level as a result of continued inbreeding. Above 
that level, as it appears to be the case on Aire and Colom, 
there is no relationship between FA and homozygosis. In 
any case, more research must be done in this respect since 
our results are based on just four insular populations. 

Besides the effect of genetic stress, insular populations 
must also withstand great environmental stress when they 
CO -inhabit their islets with the mentioned competitor- 
predator species. The gull, Larus michahellis (Naumann, 
1840), is present on Aire and Colom (Tab. 1), and the liz¬ 
ards from these two localities presented the lowest CFA2 
values. Presumably, gulls are not responsible for differenc¬ 
es in the development between islands as the GLM model 
suggests. Our results are in line with some recent papers 
that argued that their negative impact on the population 
trends of most species is not so straightforward (e.g., Oro 
& Marti'nez-Abrai'n 2007). 

The ship rat, Rattus rattus (Linnaeus, 1758), is an alien 
species that could alter endemic populations by predating 
upon them, but also because they compete for resources 


(Case & Bolger 1991). Tosqueta and Binicodrell showed 
different levels of FA, and the most major ecological dif¬ 
ference between them is the presence of R. rattus on the 
former islet (Tab. 1). Reptiles exhibit a high vulnerability to 
rodent presence (Towns et al. 2006). Our results obtained 
from the four studied populations support those of Kris- 
TENSEN et al. (2006). Thus, Tosqueta individuals could re¬ 
spond more sensitively to the presence of rats than Colom 
ones do, probably because Colom lizards have a higher lev¬ 
el of heterozygosis (Perez-Mellado et al. 2008). 

In summary, although genetic variability is considered 
the best predictor of the degree of asymmetry in popula¬ 
tions, other stressors may modulate this effect, and in our 
case it is perhaps the presence of rats. If parasites are in¬ 
deed a cause of asymmetry, it is possible that the average 
intensity of parasite load in the populations investigated 
by us does not affect the development of individuals. In 
summary, the only conclusion is that FA and parasitism 
are not related at population level in P. lilfordi on the islets 
examined. Palmer (2002) regards FA as a poor indicator 
of the quality of individual development, as there may be 
random deviations from symmetry. According to Ward et 
al. (1998), perhaps the relationship between parasitism and 
FA is not as general as suggested by other authors (Moller 
1996, PoLAK 1997, Moller 2006). There has been a recent 
debate on the importance of publication biases in our un¬ 
derstanding of ecological processes (Moller & Jennions 
2001 and references therein). Fanelli (2012a, b) argued 
that papers with negative results were disappearing from 
most disciplines in all animal and plant sciences. This is 
even more serious when dealing with relatively recent is¬ 
sues Such as FA, where statistical methodology is being 
updated continuously. Ward et al. (1998) warn “that any 
interpretation based solely on the proportion of studies re¬ 
porting a positive association maybe misleading” and can¬ 
not but agree more. 

Acknowledgements 

We thank our lab partners Ana Perez-Cembranos, Abraham 
MencIa, and Zaida Ortega, who contributed with their knowl¬ 
edge. Isabel Roldan patiently helped us year after year. Lizards 
were captured under special permits from the Servei de Proteccio 
d’Especies, ConseUeria de Medi Ambient, Balearic Government 
(permit 2011), and data were collected in accordance with its animal 
welfare standards and protocols. Financial support was provided to 
Mario Garrido by a Junta de Castilla y Leon predoctoral grant 
(Order EDU/1878/2006). This work was supported by the research 
projects CGL2009-12926-C02-02 and CGL2012-39850-C02-02 
from the Spanish Ministry of Science and Innovation. 

References 

Allendorf, E W. & G. Luikart (2007): Conservation and the 
Genetics of Populations. - Blackwell, Malden. 

Bancila, R., I. VAN Gelder, E. Rotteveel, J. Loman & J. W Arn- 
TZEN (2010): Fluctuating asymmetry is a function of population 
isolation in island lizards. - Journal of Zoology, 282(4): 266-275. 


152 


Factors determining fluctuating asymmetry in an insular lizard 


Bergstrom, C. A. & T. E. Reimchen (2005): Habitat-dependent 
associations between parasitism and fluctuating asymmetry 
among endemic stickleback populations. - Journal of Evolu¬ 
tionary Biology, 18: 938-948. 

Biaggini, M., F. Berti & C. Corti (2009): Different habitats, dif¬ 
ferent pressures? Analysis of escape behaviour and ectopara¬ 
site load in Podarcis sicula (Lacertidae) populations in differ¬ 
ent agricultural habitats. - Amphibia-Reptilia, 30: 453-461. 

Brown, C. R. & M. B. Brown (2002): Ectoparasites cause in¬ 
creased bilateral asymmetry of naturally selected traits in a 
colonial bird. - Journal of Evolutionary Biology, 15(6): 1067- 
1075. 

Case, T. J. & D. T. Bolger (1991): The role of introduced species 
in shaping the distribution and abundance of island reptiles. - 
Evolutionary Ecology, 5(3): 272-290. 

Crawford, D. J., E. Ruiz, T. F. Stuessy, E. Tepe, R Aqeveque, 
F. Gonzalez, R. J. Jensen, G. J. Anderson, G. Bernardello, 
C. M. Baeza, U. Swenson & M. Silva (2001): Allozyme diver¬ 
sity in the endemic flowering plant species of the Juan Fernan¬ 
dez Archipelago, Chile: ecological and historical factors with 
implications for conservation. - American Journal of Botany, 
88: 2195-2203. 

Crawley, M. J. (2007): The R book, 1" Edition. - John Wiley and 
Sons Ltd, Chichester, UK. 

Crnobrnia-Isailovic, j., I. Aleksic & D. Beiakovic (2005): 
Fluctuating asymmetry in Podarcis muralis populations from 
Southern Montenegro: detection of environmental stress in 
insular populations. - Amphibia-Reptilia, 26:149-158. 

Crnokrak, R & D. A. Roff (1999): Inbreeding depression in the 
wild. - Heredity, 83: 260-270. 

Cohen, J. (1988): Statistical power analysis for the behavioral sci¬ 
ences, 2““* Edition. - Academic Rress, NY. 

Fanelli, D. (2012a): Rositive results receive more citations, but 
only in some disciplines. - Scientometrics, 94(2): 701-709. 

Fanelli, D. (2012b): Negative results are disappearing from most 
disciplines and countries. - Scientometrics, 90(3): 891-904. 

Frankham, R. (1998): Inbreeding and extinction: island popula¬ 
tions. - Conservation Biology, 12: 665-675. 

Frankham, R. (2005): Genetics and extinction. - Biological 
Conservation, 126:131-140. 

Garrido, M. & V. Rerez-Mellado (2013a): Rrevalence and in¬ 
tensity of blood parasites in insular lizards. - Zoologischer 
Anzeiger, 252(4): 588-592. 

Garrido, M. & V. Rerez-Mellado (2013): Ratterns of parasitism 
in insular lizards: Effects of body size, condition and resource 
availability. - Zoology, 116(2): 106-112. 

Grafen, a. & R. Hails (2002): Modern statistics for the life sci¬ 
ences, i”* Edition. - Oxford University Rress, Oxford. 

Hothorn, T., F. Bretz & E Westfall (2009): Multcomp: Simul¬ 
taneous Inference for General Linear Hypotheses. - http:// 
CRAN.R-project.org/package=multcomp, R package version 
1.0-7. 

Ji, X., Q. B. Qiu & C. H. Diong (2002): Influence of incubation 
temperature on hatching success, energy expenditure for em¬ 
bryonic development, and size and morphology of hatchlings 
in the oriental garden lizard, Calotes versicolor (Agamidae). - 
Journal of Experimental Zoology, 292(7): 649-659. 


KrISTENSEN, T. N., R. S 0 RENSEN, K. S. Redersen, M. KrUH 0 FFER 
& V. Loeschcke (2006): Inbreeding by environmental inter¬ 
actions affect gene expression in Drosophila melanogaster. - 
Genetics, 173:1329-1336. 

Kruuk, L. E. B., j. Slate, J. M. Remberton & T. H. Clutton- 
Brock (2003): Fluctuating asymmetry in a secondary sexual 
trait: no associations with individual fitness, environmental 
stress or inbreeding, and no heritability. - Journal of Evolu¬ 
tionary Biology, 16: 101-113. 

Lailvaux, S. F. & D. J. IRSCHICK (2006): No Evidence for Female 
Association with High-Rerformance males in the Green Anole 
Lizard, Anolis carolinensis. - Ethology, 112(7): 707-715. 

Leary, R. F. & F. W. Allendorf (1989): Fluctuating asymmetry as 
an indicator of stress: implications for conservation biology. - 
Trends in Ecology & Evolution, 4: 214-217. 

Leary, R.F., F. W. Allendorf, K. L. Knudsen & G. H. Thor- 
GAARD (1985): Heterozygosity and developmental stability 
in gynogenetic diploid and triploid rainbow trout. - Hered¬ 
ity, 54(21): 9-225. 

Leung, B., M. R. Forbes & D. Houle (2000): Fluctuating asym¬ 
metry as a bioindicator of stress: comparing efficacy of ana¬ 
lyses involving multiple traits. - American Naturalist, 155(1): 
101-115. 

Lopez, F. & J. MartIn (2002): Locomotor capacity and domi¬ 
nance in male lizards Lacerta monticola: a trade-off between 
survival and reproductive success? - Biological Journal of the 
Linnean Society, 77(2): 201-209. 

Losos, J. B. (2009): Lizards in an Evolutionary Tree: Ecology and 
Adaptive Radiation of Anoles. - University of Galifornia Rress, 
Harvard. 

Mac Arthur, R. H. & E. O. Wilson (1967): The theory of is¬ 
land biogeography, i“* Edition. - Princeton University Press, 
Princeton. 

Madsen, T., B. Stille & R. Shine (1996): Inbreeding depression 
in an isolated population of adders [Vipera hems). - Biological 
Conservation, 75:113-118. 

Madsen, T., M. Olsson, H. Wittzell, B. Stille, A. Gullberg, 
R. Shine, S. Andersson & H. Tegelstrom (2000): Popula¬ 
tion size and genetic diversity in sand lizards {Lacerta agilis) 
and adders {Vipera berus). - Biological Conservation, 94:257- 
262. 

MartIn, j. & P. Lopez (2000): Chemoreception, symmetry and 
mate choice in lizards. - Proceedings of the Royal Society B: 
Biological Sciences, 267(1450): 1265-1269. 

MartIn, j. & P. Lopez (2001): Hindlimb asymmetry reduces 
escape performance in the lizard Psammodromus algirus. - 
Physiological and Biochemical Zoology, 74: 619-624. 

MartIn, j. & P. Lopez (2006): Links between male quality, male 
chemical signals, and female mate choice in Iberian rock liz¬ 
ards. - Functional Ecology, 20(6): 1087-1096. 

M0LLER, A. P. (1996): Parasitism and developmental instability of 
hosts: a review. - Oikos, 77:189-196. 

M0LLER, A. P. (2006): A review of developmental instability, 
parasitism and disease infection, genetics and evolution. - In¬ 
fection, Genetics and Evolution, 6:133-140. 

M0LLER, A. P. & M. D. Jennions (2001): Testing and adjusting 
for publication bias. - Trends in Ecology & Evolution, 16(10): 
580-586. 


153 


Mario Garrido & Valentin Perez-Mellado 


M0LLER, A. P. & J. T. Manning (2003): Growth and developmen¬ 
tal instability. - The Veterinary Journal, 166:19-27. 

M0LLER, A. P. & J. P. Swaddle (1997): Asymmetry, Developmen¬ 
tal Stability and Evolution, 1" Edition. Oxford Series in Ecol¬ 
ogy and Evolution. - Oxford University Press, Oxford. 

Oro, D. & A. MartInez-AbraIn (2007): Deconstructing myths 
on large gulls and their impact on threatened sympatric water- 
birds. - Animal Conservation, 10:117-126. 

Palmer, A. R. (1994): Fluctuating asymmetry: a primer. - pp. 
335-364 in: Marrow, T., (ed.): Developmental instability: its 
origins and evolutionary implications. - Kluwer Academic 
Publishers, Dordrecht. 

Palmer, A. R. (1999): Detecting publication Bias in Meta-ana- 
lyses: A Case Study of Fluctuating Asymmetry and Sexual Se¬ 
lection. - American Naturalist, 154: 220-233. 

Palmer, A. R. (2002): Chimpanzee right-handedness reconsid¬ 
ered: Evaluating the evidence with funnel plots. - American 
Journal of Physical Anthropology, 118(2): 191-199. 

Palmer, A. R. & C. Strobeck (1986): Fluctuating asymmetry: 
measurement, analysis, patterns. - Annual Review of Ecology, 
Evolution, and Systematics, 17: 391-421. 

Palmer, R. A. & C. Strobeck (2003): Appendix V. Fluctuating- 
asymmetry analysis: A step-by-step example. - in: Polar, 
M. (ed.): Developmental Instability (DI): Causes and Conse¬ 
quences. - Oxford University Press, Oxford. 

Perez-Mellado, V. (2002): Eradication of rodents and feral 
cats on islands. - in: Gherardi, E, C. Corti & M. Gualtieri 
(Eds): Biodiversity Conservation and Habitat Management. 
Encyclopedia of Life Support Systems. Volume III. 16. Natural 
Resources Policy and Management. - Publishing-Eolss Pub¬ 
lishers, Oxford. 

Perez-Mellado, V, J. A. Hernandez-Estevez, T. GarcIa- 
DfEZ, B. Terrasa, M. M. Ramon, J. Castro, A. Picornell, 
J. MartIn-Vallejo & R. Brown (2008): Population density 
in Podards lilfordi (Squamata, Lacertidae), a lizard species en¬ 
demic to small islets in the Balearic Islands (Spain). - Am¬ 
phibia-Reptilia, 29: 49-60. 

PiANKA, E. R. (1986): Ecology and Natural History of Desert Liz¬ 
ards. - Princeton University Press, Princeton. 

Polar, M. (1997): Parasites, fluctuating asymmetry, and sexual 
selection. - pp. 246-276 in: Beckage, N. (ed.): Parasites and 
Pathogens: Effects on Host Hormones and Behavior. - Chap¬ 
man and Hall, NY. 

Qualls, C. P. & R. M. Andrews (1999): Cold climates and the 
evolution of viviparity in reptiles: cold incubation tempera¬ 
tures produce poor-quality offspring in the lizard Sceloporus 
virgatus. - Biological Journal of the Linnean Society, 67: 353- 
376. 

R Development Core Team (2010): R: A language and environ¬ 
ment for statistical computing. - R Foundation for Statisti¬ 
cal Computing, Vienna, Austria. ISBN 3-900051-07-0, URL 
http://www.R-project.org/. 

Sarre, S. (1996): Habitat fragmentation promotes fluctuating 
asymmetry but not morphological divergence in two geckos. 
- Researches on Population Ecology, 38: 57-64. 

Sarre, S. & J. D. Dearn (1991): Morphological variation and fluc¬ 
tuating asymmetry among insular populations of the sleepy 
lizard, Trachydosaurus rugosus Gray (Squamata, Scincidae). - 
Australian Journal of Zoology, 39: 91-104. 


SoKAL, R. R. & E J. Rohlf (1995): Biometry: the principles and 
practice of statistics in biological research. - WH Freeman, 
NY. 

Soule, M. E. (1967): Phenetics of natural populations. 11 . Asym¬ 
metry and evolution in a lizard. - American Naturalist, 101: 
141-160. 

Tabachnick, B. G. & L. S. Fidell (2000): Using Multivariate Sta¬ 
tistics, 4* Edition. - Allyn & Bacon, Needham Heights. 

Towns, D. R., I. A. E. Atkinson & C. H. Daugherty (2006): 
Have the harmful effects of introduced rats on islands been 
exaggerated? - Biological Invasions, 8: 863-891. 

VAN Dongen, S. & L. Lens (2000a): The evolutionary potential 
of developmental instability. - Journal of Evolutionary Biolo¬ 
gy, 13: 326-335. 

VAN Dongen, S. & L. Lens (2000b): Correspondence: Symmetry, 
size and stress. - Trends in Ecology & Evolution, 15: 330-331. 

VAN Valen, L. (1962): A study of fluctuating asymmetry. - Evolu¬ 
tion, 16:125-142. 

Vervust, B., S. van Dongen, I. Grbac & R. van Damme (2008): 
Fluctuating asymmetry, physiological performances, and 
stress in island populations of the Italian wall lizard [Podards 
sicula). - Journal of Herpetology, 42: 369-377. 

Ward, D. E, F. Thomas & R. Poulin (1998): Fluctuating asym¬ 
metry and parasitism in six New Zealand insects. - Acta 
Oecologica, 19(5): 409-412. 

Warner, A. & R. Shine (2006): Morphological variation does not 
influence locomotor performance within a cohort of hatchling 
lizards (Amphibolurus muricatus, Agamidae). - Oikos, 114(1): 
126-134. 

Watson, P. J. & R. Thornhill (1994): Fluctuating asymmetry and 
sexual selection. - Trends in Ecology & Evolution, 9: 21-25. 

Wiener, J. G. & P. J. Rago (1987): A test of fluctuating asymmetry 
in Bluegills (Lepomis macrochirus Rafinesque) as a measure of 
pH-related stress. - Environmental Pollution, 44(1): 27-36. 

Young, A., T. Boyle & T. Brown (1996): The population genetic 
consequences of habitat fragmentation for plants. - Trends in 
Ecology & Evolution, 11: 413-418. 

Zar, j. H. (2010): Biostatistical analysis, 5* Edition. - Prentice 
Hall, Upper Saddle River. 

Zhdanova, N. P. & V. M. Zakharov (2006): Developmental sta¬ 
bility of sand lizard [Lacerta agilis L.) under optimal and non- 
optimal conditions. - pp. 209-220 in: Corti, C., P. Lo Cascio 
& M. Biaggini (eds): Mainland and insular lizards. A Medi¬ 
terranean perspective. - Firenze University Press, Firenze. 

ZuuR, A. F, E. N. Ieno & C. S. Elphick (2010). A protocol for 
data exploration to avoid common statistical problems. - 
Methods in Ecology and Evolution, 1(1): 3-14. 


154 


